Efimov physics in heteronuclear four-body systems 
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We study three- and four-body Efimov physics in a heteronuclear atomic system with three 
identical heavy bosonic atoms and one light atom. We show that exchange of the light atom 
between the heavy atoms leads to both three- and four-body features in the low-energy inelastic 
rate constants that trace to the Efimov effect. Further, the effective interaction generated by this 
exchange can provide an additional mechanism for control in ultracold experiments. Finally, we find 
that there is no true four-body Efimov effect — that is, no infinite number of four-body states in 
the absence of two- and three-body bound states — resolving a decades-long controversy. 

PACS numbers: 



Few-body physics has benefitted greatly from ultracold 
experiments on quantum gases in recent years. Many 
long-standing predictions have been verified and new 
phenomena identified due to the extraordinary ability to 
control and measure these systems. One of the more 
spectacular results for few-body physics was the obser- 
vation of Efimov physics in Cs [1], experimentally vali- 
dating Efimov's surprising, 35-year old prediction [2] by 
confirming its predicted influence on three-body recom- 
bination [3]. Several measurements of Efimov physics 
have now been made [1, 4-13] and even include features 
traceable via theory [14] to four-body processes [10, 15]. 

One natural question to ask is whether there is an Efi- 
mov effect for A>3 bodies. The Efimov effect in this 
case is defined as the existence of an infinite number of 
A-body bound states when no subsystems are bound. 
Part of the answer was provided in Ref. [16]: there is 
no Efimov effect for A>3 equal mass particles. While 
generally interpreted to imply no Efimov effect is possi- 
ble for A>3 (see, for example, [17]), this result does not 
preclude the possibility of an Efimov effect for systems 
with unequal masses. In fact, Ref. [18] recently reported 
that three identical fermions interacting resonantly with 
a fourth particle do have an Efimov effect for a small 
range of mass ratios. The question of an Efimov effect 
for A=4 with three identical bosons, however, remains 
open despite prior study: H 3 L systems (H and L arc 
heavy and light particles, respectively) were analyzed in 
Ref. [19] with the conclusion that no Efimov effect occurs. 
When this system was revisited in Ref. [20], however, the 
opposite conclusion was reached. 

In this Letter, we settle this controversy: there is no 
true four-body Efimov effect for H 3 L with bosonic Hs. 
Although our conclusion is in agreement with Ref. [19], 
our reasoning is very different. We have, however, identi- 
fied one universal four-body state attached to each three- 
body Efimov threshold. We also show that the low- 
energy scattering observables simultaneously display dis- 
tinct three- and four-body features characteristic of the 
three-body Efimov effect. Moreover, we find that, in the 
context of ultracold collisions, the s-wave two-body scat- 



tering length anL between H and L atoms can be used 
to tune the effective heavy-heavy scattering length a* HH 
when auL > 0, opening up new avenues for control in 
few-body systems. 

Since a direct solution of the four-body problem 
remains a substantial challenge, we apply the Born- 
Oppcnhcimer (BO) approximation as described in 
Refs. [19, 20], assuming that the mass m of L is much 
less than the mass M of H . Although we will focus on 
bosonic Hs, our analysis can be generalized straightfor- 
wardly. 

We first recall that the three-body Efimov effect oc- 
curs when at least two scattering lengths are much larger 
than a characteristic two-body interaction range tq [2]. 
Thus, for H2L, Efimov physics will occur for laHL^fo- 
Applying the BO approximation with p the H+H dis- 
tance yields a BO potential with the universal long- 
range behavior — xffi 2 /2mp 2 , x ~0. 567143, in the region 
ro<Cp^;|a£fL| [17, 21]. The three-body problem has thus 
been reduced to an effective two-body problem whose 
interaction is the Efimov potential, giving the character- 
istic geometric relation between bound state energies 
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where Sq=XqM /2m — 1/4. An infinite number of Efimov 
states are thus obtained in the limit |<Zhz,|— >oo, corre- 
sponding to a zero-energy HL bound state. 

Similarly, all of the universal results [17] for low-energy 
three-body scattering apply to H^L. For instance, when 
Ohl>0, there is a weakly bound HL molecule, and 



a* HH cx cot[s ln(oHL/Vo) + $]a>HL 



(2) 



where r$ is the characteristic range of the two-body in- 
teraction. The poles in a* HH occur at those values of auL 
when an Efimov state becomes bound, so their positions 
are characteristically log-periodic in a^L- Their over- 
all position, however, is determined by the short-range 
three-body phase $ [3, 22]. Equation (2) thus shows that 
aHL provides a control over a* HH that could prove advan- 
tageous when a* HH cannot easily be controlled directly 
via, for example, a Feshbach resonance [23]. 
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In general, the Hs also interact directly. However, so 
long as the direct H+H interaction is short-ranged - 
even if it is repulsive — the large-p behavior of the ef- 
fective H+H interaction is unchanged, and the Efimov 
effect remains. A direct H+H interaction changes only 
Eq (or, equivalently, $) and not the characteristic log- 
periodic behavior. Consequently, manipulating the di- 
rect H+H interaction allows control of $ and thus the 
position of the family of log-periodic Efimov features. 

Applying the BO approximation to H 3 L reduces it to 
an effective three-body problem for H 3 . We thus expect 
that our knowledge of Efimov physics and universality for 
three identical bosons should apply and provide at least 
a basic understanding of the system. For simplicity, we 
assume that the H atoms do not interact directly. Per 
the argument above, we know this will not materially 
affect our conclusions. The H+H interaction then comes 
solely from mediation by the L and is characterized by 
the effective H+H scattering length a* HH from Eq. (2). 
Based on the result for three bosons, we expect an Efimov 
effect for H 3 when \a* HH \^oo [20]. 

To determine whether the \a^ H \-^oo limit actually 
produces an Efimov effect for H 3 L, we adopt the fol- 
lowing definition [16, 19]: a true four-body Efimov ef- 
fect exists if the four-body system possesses an infin- 
ity of stable bound states when there is a zero-energy 
three-body bound state and no other two- or three-body 
bound states. Under these conditions, the four-body Efi- 
mov effect is completely independent of any three-body 
Efimov effect, making the resulting states fundamentally 
different from the universal four-body states discussed in 
Refs. [14, 24-26]. 

If it exists, a true four-body H 3 L Efimov effect will 
occur for la^]— »-oo and clhl<0- This case gives a zero- 
energy H2L bound state but no HL bound state and 
is precisely the case identified in Ref. [19]. They con- 
cluded that no Efimov effect is possible because the BO 
potential surface has no long-range component. Naus 
and Tjon [20] correctly pointed out that it is not the BO 
potential surface that must behave as — i?~ 2 , but rather 
the adiabatic hyperspherical potential where the hyper- 
radius R measures the overall size of the system (see, for 
example, Refs. [27, 28] for a discussion of hyperspherical 
coordinates in this context). They concluded that be- 
cause \a* HH \~^oo : the conditions of the Efimov effect for 
H 3 are fulfilled and there is thus a four-body Efimov ef- 
fect. Unfortunately, neither analysis provides a definitive 
answer. 

Although otherwise sound, the argument of Ref. [20] 
is limited by its reliance on the BO approximation — a 
problem they identified but did not address. The issue 
with the BO approximation is that for uhl<0, L becomes 
unbound when p> | ajyz, | [19], and its spectrum becomes 
continuous. Since a zero-energy L is no longer fast com- 
pared to the Hs, the BO approximation breaks down. 
Naus and Tjon followed Ref. [19] and simply set the BO 
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FIG. 1: Schematic energy spectrum for H3L. Solid lines 
denote bound states and dashed lines denote resonances. 



potential at such points to zero. Under this questionable 
approximation, a* HH can still be defined and adjusted to 
infinity, so that they conclude there is an infinite series 
of bound H 3 L states, i.e. the four-body Efimov effect. 

Where the BO approximation fails, however, the adia- 
batic hyperspherical approximation is perfectly valid. In 
particular, the lowest H2L adiabatic hyperspherical po- 
tential, which corresponds to the BO potential, crosses 
zero energy at a hypcrradius R^\o,hl\ much like the BO 
potential, but remains discrete and increases to a bar- 
rier with height ~~i-/ajj L [29]. For R beyond this barrier, 
the potential falls to zero as l e g(l e ff + l)h 2 /MR 2 with 
U=3/2 [30]. 

Having taken care to define the H2L hypcrradius R 
such that it reduces to the H+H distance p in the 
M/m— >oo limit, we can regard this adiabatic hyper- 
spherical potential as an effective H+H interaction. 
Since the Efimov effect for the three Hs can only oc- 
cur for s-wave two-body interactions, i.e. 4ff=0, the fact 
that Z e ff=3/2 prohibits the Efimov effect. This conclu- 
sion has been confirmed by direct calculations of the 
four-body adiabatic hyperspherical potentials for H 3 L 
with M/m=30 [31] using the correlated Gaussian ap- 
proach [32]. 

Even though our finding no true four-body Efimov ef- 
fect agrees with the conclusions of Ref. [19] and contra- 
dicts the conclusions of Ref. [20], we believe the former 
were right for the wrong reason and the latter underes- 
timated the consequences of the breakdown of the BO 
approximation. In the end, it is this breakdown that ex- 
cludes the possibility of a true four-body Efimov effect in 
this system. 

The possibility of true four-body Efimov states is not 
the only phenomenon of interest in H 3 L. Continuing to 
o-hl>0 such that aHL^^o, the BO approximation dis- 
plays no pathologies since L is bound for all H config- 
urations, and we can safely think about the three-body 
H 3 motion on the lowest BO potential surface. Based on 
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the known three-body results [17, 27, 33], we expect an 
adiabatic hypcrspherical potential of the form 
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with <7o=l-00624 and /j,=M/y/3 that approaches the 
three-body break-up threshold. In this case, that thresh- 
old corresponds physically to HL+H+H. The lower limit 
of R in Eq. (3) is modified from the usual three-body 
problem [33] due to the fact that the characteristic range 
of the effective two-body H+H potential is no longer ro, 
but rather auL as defined by the size of the HL bound 
state. 

When anL is tuned to give a HH —>oo, the Efimov poten- 
tial (3) extends to infinity, producing an infinite series of 
four-body bound states below the HL+H+H threshold 
with binding energies E n+1 /E n =e- 27 <l 9a . These states 
are not true four-body Efimov states since there is an 
HL bound state, but they can be regarded as three-body 
Efimov states of HL+H+H . They arc indicated with 
the notation (HL)H 2 in Fig. 1 where the H 2 L+H and 
HL+H+H thresholds intersect since a*^— >oo at these 
points. 

Figure 1 sketches the energy trajectories for H3L as 
a function of anL- In addition to the (HL)H 2 Efimov 
states, there are the H 2 L Efimov states, and associated 
with each of these we find — for both M/ra=50 and 
M/m=30 — one H3L state that appears to be the analog 
of the universal four-boson states in Refs. [14, 15, 24-26]. 
We find that its binding energy is universally related to 
the binding energy of the associated H 2 L Efimov state 
by Eh 3 l/ Eh 2 l~QA for both of the mass ratios we have 
calculated. Note that the binding energies are defined 
relative to the next lowest breakup threshold (HL+H+H 
for H 2 L states and H 2 L+H for H3L states). 

Figure 2 shows our numerically calculated three-body 
recombination rates K3 for HL+H+H^rH 2 L+H with 
M/m=30 (see Ref. [28] for details of our numerical meth- 
ods). Although there may be several final H 2 L Efimov 
states available, our calculation shows that recombina- 
tion into the most weakly bound Efimov state dominates. 
In fact, the main peaks in Fig. 2(c) occur where an H 2 L 
Efimov state just becomes bound. The separation be- 
tween the main peaks is thus determined by sq from the 
H 2 L Efimov effect: for M/m=30, a^/a^=e 7r / s °=4.34. 
Note that the adiabatic hypcrspherical approach gives 
e f/«o = 3.96 [33, 34], giving an indication of the BO ap- 
proximation error for this mass ratio. For consistency, we 
will quote only BO results in the rest of this Letter. The 
factor of 5.5 between these two main peaks in Fig. 2 does 
not match this prediction because the criterion a/a3>ro 
is not well satisfied. 

Since each main peak corresponds to a pole of a* HH , 
K3 shows Efimov features characteristic of the H3 mo- 
tion where \a* HH \^>aHL- When this condition is satisfied, 
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FIG. 2: K 3 for HL+H+H^H 2 L+H when M/m=30. The 
main peaks in (c) reflect the Efimov physics of H2L and occur 
where (d) a* HH diverges. The expanded plots (a) and (b) high- 
light the features on each main peak related to H3L states. 
The vertical dashed lines indicate where a* HH =oo. 



the rates are given by the usual universal three-body ex- 
pressions [17, 35] with the short-range length scale set to 
o-hl 



C 



sinh(2?y) \a 



HH I 



K 



(a' HH >0) 



V sin 2 [g ln(\a* HH \/a HL ) + <l> / }+ sinh 2 (77) 
C 

- sm 2 [g ln(a* HH /a HL ) + $] (a* HH ) 4 . (4) 



In these expressions, C and C are universal constants. 
But, because the final H 2 L state is an Efimov state, $, $' 
and r\ depend not on short-range four-body physics, but 
rather on the short-range physics of the H 2 L states - 
no additional four-body parameter is needed [14, 24, 25]. 
The a HH -K 3 projection in Fig. 3 supports this conclu- 
sion, showing that the positions of the Efimov features 
related to the H3 motion described by Eq. (4) are ap- 
proaching a universal position in the limit Uhl^to- 

Interestingly, for larger mass ratios, H 2 L Efimov states 
with non-zero orbital angular momentum j are possi- 
ble [33, 34]. In this case, the universal constant sq in 
Eq. (1) is determined from Sq=-)^M /2m— + — 1/4, 
and there will be an Efimov effect for H 2 L so long as 
Sq>0. Higher angular momentum Efimov states have 
not yet been observed because experiments have focused 
mostly on identical particles which have no such states 
and because j>0 Efimov states produce extremely nar- 
row recombination peaks as seen in Fig. 3. The mass 
ratio in Fig. 3, M/m=50, supports j=2 Efimov states, 
but is not sufficient for higher j states. Because Z e g —j^0 
for the H 2 L system, there is no (HL)H 2 Efimov effect as- 
sociated with these states. Their narrow K3 peaks thus 
show no substructure of the sort seen on the main peaks. 

The periods for each of the three different families of 
Efimov peaks are indicated in Fig. 3. Because none of the 



FIG. 3: K 3 for HL+H+H^H 2 L+H with M/m=50 simul- 
taneously showing the dependence on auL and on a* HH . The 
H 2 L Efimov physics is highlighted by the cihl-Ks projection 
while the a* HH -Kz projection illustrates the (HL)H 2 Efimov 
physics. Note that a* HH increases into the figure. 



respective scattering lengths are strongly in the univer- 
sal limit, however, the calculated periods do not match 
the predicted ones. The main j=0 H2L Efimov peaks 
highlighted by the CLHL-K3 projection should have a pe- 
riod of e n / s °=3.08. The calculated spacings are larger 
than this, but appear to be approaching the expected 
value as clhl increases. Similarly, for j=2, the expected 
period is 10.46 while the calculated one is 17.5, and the 
predicted period of the (HL)H 2 Efimov substructure on 
the main peaks highlighted in the a* HH -K^ projection is 
e*'/ 90 =22.7. The magnitude of the deviations from the 
predictions likely reflects what can be observed since it 
is difficult to penetrate deeply into the universal regime 
experimentally. 

So far, we have not taken advantage of all of the free- 
dom that this hetcronuclear system affords to manipulate 
the Efimov features. In particular, since a* HH depends 
only on the total H+H interaction, i.e. effective plus di- 
rect interactions, it can be tuned by either interaction — 
or both. Intriguingly, this freedom also allows experimen- 
tal control of both a* HH and $ (or $'). Controlling the 
former via either clhl or a direct interaction with scatter- 
ing length cihh allows the various Efimov features seen, 
for instance, in Figs. 2 and 3 to be mapped out. Con- 
trolling the latter makes it possible to shift all of these 
features — something not possible so far in three-body 
systems. Moreover, in the neighborhood of a pole in a* HHl 
it should be possible to exert both types of control largely 
independently. 

To illustrate the effect of tuning the direct interaction, 
we show in Fig. 4 the numerically calculated rates for 
the relaxation process H2L(n)+H^H2L(n — 1)+H as a 
function of clhh where n labels the most weakly bound 
H2L state. This tuning was accomplished by including 
a short-range, direct H+H interaction in addition to an 



FIG. 4: V ie i for H 2 L{n)+H^H 2 L{n - 1)+H for (a) n=l 
and (b) n=2. The dashed lines indicate the positions where 



effective H+H interaction. The behavior of V le \ in Fig. 4 
is also found in the three-boson system. In fact, when 
a *HH ^> clhl > 0, V Te i has the same form as in three- 
boson systems [17] 



A sinh 2ry 



sin 2 [g \n(a* HH / cihl ) + + sinh 2 r\ 



l HHi 



(5) 



where A is a universal constant, but 77 and $ depend on 
the short-range details of the relaxed H 2 L bound state 
— although not on a separate four-body parameter. 

To conclude, we have studied Efimov physics in the 
four-body hetcronuclear system H3L with bosonic H 
atoms, showing that there is no true four-body Efimov 
effect. We have, however, identified a universal four- 
body state linked to each H 2 L Efimov state. Perhaps 
more significantly, we have shown that the richness of 
the heteronuclear system can be exploited to reveal a va- 
riety of Efimov features with different universal scaling 
by scanning just am,. That same richness allows more 
opportunities for controlling the system as well, which 
is potentially very important for experiments seeking to 
study or utilize these processes. For instance, taking ad- 
vantage of the BO approach introduced in Ref. [19], we 
pointed out that scanning clhl gives a means for control- 
ling the scattering length between the heavy atoms — a 
tool that could prove useful experimentally. Heteronu- 
clear four-body and larger systems thus provide consid- 
erable possibilities for studying Efimov physics that were 
not previously anticipated. 
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